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The Weyl semimetal is characterized by three-dimensional linear band touching points called Weyl 
nodes. These nodes come in pairs with opposite chiralities. We show that the coupling of circularly 
polarized photons with these chiral electrons generates a Hall conductivity without any applied 
magnetic field in the plane orthogonal to the light propagation. This phenomenon comes about 
because with all three Pauli matrices exhausted to form the three-dimensional linear dispersion, the 
Weyl nodes cannot be gapped. Rather, the net influence of chiral photons is to shift the positions 
of the Weyl nodes. Interestingly, the momentum shift is tightly correlated with the chirality of the 
node to produce a net anomalous Hall signal. Application of our proposal to the recently discovered 
TaAs family of Weyl semimetals leads to an order-of-magnitude estimate of the photoinduced Hall 
conductivity which is within the experimentally accessible range. 

PACS numbers: 73.43.-f, 03.65.Vf, 72.40.+w 


Introduction - Previously thought to be an asset ex¬ 
clusive to massless elementary particles such as photons, 
chirality has now become a defining emergent property 
of electrons in such crystalline materials like graphene [l| 
and the surface of topological insulators [H. An excit¬ 
ing new addition to the growing list of “chiral electronic 
materials” is the Weyl semimetal [il3- In this three- 
dimensional analogue of graphene, the Weyl nodes act 
as magnetic monopoles in momentum space [H . In sharp 
contrast to the two Dirac nodes in graphene that behave 
more or less as separate low-energy degrees of freedom, 
a pair of Weyl nodes carrying opposite chiralities are in¬ 
herently tied together in a non-local manner through the 
“chiral anomaly” mechanism which allows the dissipa¬ 
tionless transfer of charge between them The subject 
has undergone vigorous research lately due to reports of 
material predictions, followed immediately by their syn¬ 
thesis and confirmation of their Weyl characters by pho¬ 
toemission experiments and non-local transport measure¬ 
ments. A striking example of success along this line in 
recent years is the transition metal monophos phi de fam¬ 
ily that includes TaAs, TaP, NbAs, and NbP lOmij. 


In light of the rapid maturity of the Weyl semimetal 
research, one should ask whether the unusual electronic 
transport properties in both the static and dynamic 
regimes of the Weyl semimetal can be exploited in high 
speed electronics or to provide a new means of reveal¬ 
ing dynamic topological effects It was recently 

suggested that exposing a two-dimensional Dirac mate¬ 
rial to a circularly polarized (CP) light could lead to a 
novel form of Hall effect due to the effective gap opening 
at the Dirac point 21, 22|. A natural question arises as 
to what happens when the three-dimensional, linearly 
dispersing band of electrons couple to an intense CP 
light, as schematically presented in Fig. [H and whether 


a non-trivial Hall effect can be induced. Here we present 
a direct consequence of coupling CP light to the Weyl 
fermion: the appearance of the anomalous Hall effect 
(AHE). 

In this work, we show that even when the AHE 
is absent due to symmetry reasons in typical Weyl 
materials, it must be generally present in all Weyl 
systems when coupled to a CP light source. Based on a 
low-energy effective Weyl-Eloquet Hamiltonian analysis 
and symmetry considerations, we discuss the induced 
AHE as a generic and readily observable phenomenon 
in Weyl semimetals. We further demonstrate this effect 
using a concrete microscopic model, and apply our study 



FIG. 1. (Color online) Schematic figure for a driven Weyl 
semimetal. Blue and red circles indicate Weyl nodes with 
opposite chiralities Xw' "Fhe node positions are shifted by 
the chiral photons in a chirality-dependent manner and the 
shift is proportional to the square of the driving amplitude A. 
As a result, even though the total momentum shift along the 
driving direction is zero Aki^^ = 0), the overall Chern 
vector shift has a finite ^-component ^ 

0), resulting in a photoinduced anomalous Hall conductivity 
axy (see Eq. (|6|)). 
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to the TaAs family of Weyl semimetals with an estima¬ 
tion of the experimental feasibility. Other interesting 
works have addressed the driven three-dimensional Dirac 
electrons ( 2 ^. [ 2 ^ or gapped systems [IBI, [ 2 ^ . However, 
neither systems enjoy the unique topological protection 
of the Weyl nodes, whose response, under a CP beam, 
results in the experimental signature of the AHE as we 
explain below. 


General tight-binding model analysis: The low-energy 
Hamiltonian for a single Weyl node can be parameterized 
most generally by a real 3-vector a and a 3 x 3 real matrix 
/3: 


light when cos0 = Al, and vanishes for a linearly polar¬ 
ized beam. 

The CP light has two contributions to the Weyl- 
Floquet Hamiltonian in Eq. O- The first (/3^) term 
originates entirely from the two-band Weyl Hamiltonian 
while the second term proportional to CxC^ — 
CyCl is due to high-energy band-mixing and must be 
kept in the construction of a general Weyl-Floquet ef¬ 
fective model. For later convenience, we re-parameterize 
-{i/2){C,Cl-CyCl) = /^oCTo + with four real num¬ 
bers K,i = {—i/A)TT[cFi{CxC^y — CyCl)]. Therefore, the 
influence of the drive to A^-order on the 2x2 Weyl band 
can be summarized as: 


Hw{^ = QiOiiCro + qiPijCJj, ( 1 ) 

with a canonical 3D linear dispersion Ew{q) = c^iqi E 
' Einstein summation is assumed in 
the momentum qi and the Pauli matrix (jj. The chi¬ 
rality of the Weyl node is given by the sign of the de¬ 
terminant of yd: xw = sgn[Det(yd)]. For a practical 
band structure, this 2x2 Weyl Hamiltonian is embed¬ 
ded within a large Wband Hamiltonian as a low-energy 
sector. In order to have the most general consideration of 
the photoinduced AHE, we therefore consider a Wband 
tight-binding Hamiltonian H(k). The expansion oi H(k) 
around a Weyl node at kw to first order in g = ^ — kw 
can be put in the block form: 


Hwf{^ = {Qi - Sqi)aiao + {qi - Sqi)/3., 




61 J, ■ ao, 


Sqi = 


e'^A^Ay cos (f) 


LO 


[^jklf^xjf^yk T ^l]{P ^lii 


Sjii = —Sq • a ■ 


e^AxAy cos (p 


-no- 


( 5 ) 


Sqi and Sja represent the shifts of Weyl momenta and 
chemical potential, respectively. 

The photoinduced Weyl node shift has an immediate 
consequence on the electronic transport. According to 
Ref. [6[ , for a Weyl system where the chemical potential 
situates at the nodes (i.e. // = 0), the leading order 

change of the anomalous Hall conductivity is governed 
by the momentum shifts of every Weyl node as: 


U^H{kw + ^U = 
= 


Hnn{q) + 0{q^), 

fHw{q) QiCi \ 

V Do + qiDiJ 


( 2 ) 


where the matrix Dq gives the {N — 2) high-energy states 
at A? = kw- Eigenvectors at A? = kw are used to con¬ 
struct the unitary matrix U whose first two columns cor¬ 
responding to the two zero-energy states. In general, the 
2 X {N — 2) matrices Ci mix the Weyl bands with high- 
energy bands linearly in q and should not be ignored. 

Without loss of generality, we consider an incident elec¬ 
tromagnetic wave in the ^-direction with the vector po¬ 
tential A{t) = (Ax cos(ujt), A^i sinjujt -h 0). Following 
standard routes [2l|, the Peierls substitution Hi[^{q) 

and averagin g o ver one drive cycle leads 
to the effective Hamiltonian |27j : 


ScTij = where ^ ■ 5q^p. (6) 


Different Weyl nodes are labeled by the superscript / and 
is the change of the Chern vector Provided that 
the momentum shifts 5q^^^ among Weyl nodes of different 

chiralities Xw cancel out, one would expect the 

AHE induced at the A^-order proportional to the inten¬ 
sity of the incident beam. Note that the photoinduced 
self-doping (dji ^ 0) of the Weyl nodes contributes to 
5(Jij at a higher order 0{6jA^) = O(A^) Q of the drive 
intensity, which is negligible. If = 0, one can show 
the ^-component of the Chern vector shift to be 




e^AxAy cos 0 

UJ 



[Cof(/3(^)),,]^ 

|det(/3(^))| 


( 7 ) 


77eff(g1 = 


H^F{q) 0{q) + 0{A^) \ 

0{q) + 0{A^) Do + 0{q) + 0{A^)) 


( 3 ) 


where the 2x2 Weyl-Floquet part reads 


ifwF® = 77 w(9) - X 


UJ 


jkPxiPyjO’k • ( 4 ) 


Higher-order terms O(g^) + 0(gA^) + 0(A^) are ignored. 
Clearly, the photoinduced effect is maximized for a CP 


where is the ij element of the cofactor ma¬ 

trix of The alternating signs of x^^^ between a node 
and an anti-node, or a monopole and an anti-monopole, 
are precisely cancelled by the same sign change in 5qi^\ 
As such, each Weyl node has a positive-definite contri¬ 
bution to 5vz- Consequently, a finite Hall conductivity is 
induced in the plane perpendicular to the incident beam 
and its sign is determined by the photon chirality through 

^xy COS (j). 

In practice, the pre-drive chemical potential /i may 
not cross the Weyl node as we discussed so far. As long 
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as /i before the drive lies in the linearly dispersive regime 
near the node, one can show that the our resultant 
6aij is unaffected [27|. In fact, using a frequency of 
excitation 0 < cj < 2// for slightly electron-doped Weyl 
node would ensure no absorption, as assumed here, and 
thus should lead to a more straightforward analysis of 
the experimental results regarding the photoinduced 
AHE. 


Symmetry consideration for multiple Weyl nodes: In 
a typical Weyl material, there are multiple sets of Weyl 
nodes related by the space group and time-reversal (TR) 
symmetries [Sl, llO|. The parameters a, /3, and Ci charac¬ 
terizing symmetry-related nodes transform accordingly. 
We now derive their transformation rules and discuss how 
the nodal shifts 5q and dji among the symmetry-related 
nodes are correlated. 

A pair of Weyl points (with / = 1,2) related by 
some space-group symmetry satisfy Rk^ = k^ ^ where 
R G 0(3) defines the particular symmetry operation. 
The parameters quantifying the Weyl nodes transform 
according to 


a. 


( 2 ) 


R 

1 hij LX j , 


Cf (8) 


As a result, the chiralities of nodes related by a rotation 
(Det[R] = +1) are the same, whereas a mirror or inver¬ 
sion operation (Det[R] = — 1) produces a sign difference. 
On the other hand, nodes related by an anti-unitary sym¬ 
metry involving TR operation satisfy 


Due to the presence of ^ = diag(l, — 1,1), the 

chiralities of these two nodes are related by —Det[R]. For 
instance, the TR operation with R = diag(—1, —1, —1) 
gives the same chirality for the nodes at kw and Rkw = 
—kw- 

Some proposed Weyl semimetals break TR symme¬ 
try but preserve inversion (I). Our symmetry analysis 
predicts and 

/^(2) = ^(1) for a pair of I-related nodes and thus, using 
Eq. (O, 

(I —related) (10) 


Since these two nodes must have opposite chiralities 
(Det[R] = —1), the product is identical for 

both of them. Hence, from Eq. these two nodes 

contribute additively to the Hall conductivity. A similar 
conclusion holds for an I symmetry breaking Weyl ma¬ 
terials for which TR is a good symmetry. In this case 
we find = —dfd\ and 

for TR-related nodes, and thus 

(5^^) = (TR—related) (11) 

Since Eq. (|6]) once again predicts a con¬ 

structive contributions from the TR-related nodes. 



FIG. 2. (Color online) (a) The 12 Weyl nodes in the first 
Brillouin zone of the undriven lattice Weyl model studied in 
Ref. [ 2 ^. (b), (c), and (d) Photoinduced Weyl node shifts 
in the presence of a CP drive along the ^-direction. Nodes 
with positive (negative) chirality are represented by red (blue) 
dots, and (I) labels the I-th Weyl node (see Eq. (O). Open 
circles denote the undriven Weyl node positions. Only pro¬ 
jections on the kz = 27r, ky = 27r, and kx = 27r planes are 
shown as we increase from 0 to 0.3. The shifts in (b) 
are magnified 10 times for a better visualization. Details of 
the model calculation can be found in [^. Parameters used: 
(t,6,A,cc;) = (1,3,1,1.1). 


Lattice model analysis: To illustrate our findings, we 
now proceed to study a concrete lattice model for Weyl 
semimetal of inversion symmetry breaking type [HI . The 
model is parameterized by the hopping (t), spin-orbit 
coupling (A) and inversion breaking potential (e) [ 2 ^ . 
The bare system possesses 12 Weyl nodes labeled as (see 
Fig.Efa)): 


= (±fco, 0 , 27r), = (0, ±fco, 27r), 

k^w+ = (±^o, 27r, 0), = (27r, ±fco, 0), 


= (0,27r,±fco), 


= (27r,0,±fco), (12) 


along with their chiralities as subscripts ± and a charac¬ 
teristic momentum /cq = 2 sin“^ [e/(4A)]. In the presence 
of a CP drive, each node undergoes a momentum shift 
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given by: 


Sq 


^1.2) _ ^^ 3 , 4 ) _ 8r]{l - 7?)A^ 


00 


Jgf.e) ^ _^^-<7.8) ^ 8r?(l + r?)A^^ 

^ 00 e 

dq_ -dq+ - ^ 


(13) 


with T] = y^l — e^/16A^. The momentum shifts clearly 
correlate with the chiralities. The photoinduced self¬ 
dopings are: 


(^4) 

00 16 

and = 0. The anomalous Hall conductivity ob¬ 

tained is (up to H^-order): 


^ _ e‘^ e‘^A‘^ f QArjX^ et^ \ 

^ v:r 


(15) 


which agrees with the expectation of the low-energy anal¬ 
ysis. 

So far all our discussion have been done within 
the effective Hamiltonian scheme valid for a small 
driving amplitude. A more stringent check free 
from this assumption is offered by solving the full 
Floquet Hamiltonian whose matrix elements are 
{n'\HF\n) = Hn-n' + noo5n,n'h with n(n') denoting the 
n(n')-th Floquet band. The topological character of 
the Weyl nodes indeed protects them against the drive, 
and only their locations are shifted continuously in a 
chirality-dependent way. Figure Efb-d) demonstrate the 
Weyl node shifts obtained from diagonalizing the full 
Floquet Hamiltonian as we increase the driving field 
strength. It is obvious that nodes with opposite chirality 
have opposite shifts in kz, leading to an imbalance of 
duz and eventually, a net anomalous Hall conductivity. 
The overall field dependence of the sum of Akz weighted 
by chiralities is shown in Fig. [3fa), which exhibits the 
anticipated increment (we have put e = 1). This 
quantity equals to the Chern vector Suz in the small A 
regime (Eq. (j6j)). On the other hand, the resultant self¬ 
doping shows the expected dependence as illustrated 
in Figure lll^b). However, it only corresponds to a higher 
order effect and does not affect the photoinduced AHE. 


Application to TaAs family of Weyl semimetals: TaAs 
was recently shown to host 24 Weyl nodes in the Bril- 
louin zone [l0|, llll. Il5l . There are two nonequivalent sets 
of Weyl points [l0| that we denote as: the P-set con¬ 
sisting of 8 nodes at (±/c^, 0) and (±/c^, 0), 

and the Q-set with 16 nodes at (±/c^, and 

{±k^,±k^,±kf). The values of k^.kyAxAyA? foi* 
various TaAs family of materials can be found in Ref. US- 






FIG. 3. (Color online) The field dependence for (a) the total 
Akz weighted by chiralities and (b) the absolute self-doping 
\Sy\ averaged over 12 Weyl nodes for the driven Weyl lattice 
model. Same parameters used as in Fig. [2l 


Applying the symmetry transformation rules derived be¬ 
fore, we find a simple rule for the node shifts for both 
Weyl sets: 

1 < / < 8 , 

Sq(QJ) = ^(QX)sqQ^ 1 < / < 16, (16) 

where Sqf^^^ oc is the same for all 8(16) symmetry- 
related nodes. The ^-direction momentum shift faithfully 
follows the chirality of each Weyl node, whereas the 
same is not true for the shifts in the xy-plane. Therefore, 
the overall photoinduced AHE for TaAs family is: 

g2 

<7xy = (8 X 5q^ -h 16 x 5q^) . (17) 


Experimental Realization: To guide future exper¬ 
imental efforts to observe these effects, we discuss 
the conditions needed and likelihood of success. In the 
supplemental we detail our calculations of the magnitude 
of the Hall voltage [l^, where we assume the use of a 
CW CO 2 laser, with a roughly constant illumination 
across the contacts, spaced ^ 100 /im apart. Using 


0, 


established values for the Eermi velocity [IJ, Il6l 
reflectance, optical [ 2 ^ and DC conductivities [l8|, we 
find a 100 nm thick film of TaAs would produce 
a Hall signal ^ 130 nV at room temperature, with a 
DC current of lA and IW laser power. We note that 
the small penetration depth requires the use of thinner 
samples to ensure most of the current is modulated by 
the light (the induced voltage is inversely proportional 
to the square of the thickness). While such signals 
should be straight forward to detect, even higher values 
may be possible through the use of pulsed lasers with 
higher peak electric fields, with the Hall signal detected 
through Earaday rotation measurements. 


Conclusion: An assortment of analysis based on Elo- 
quet theory is carried out on models of Weyl semimetals 
to argue that the AHE is induced generically by apply¬ 
ing an AC electromagnetic field of a definite chirality. 
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in the plane orthogonal to the incident beam. This pho- 
toinduced AHE originates from the remarkable monopole 
nature of Weyl nodes. The induced Hall conductivity 
scales linearly and continuously with the field intensity. 
The nodal shift itself may be observable by the pump- 
probe ARPES, while the photoinduced AHE can be de¬ 
tected via DC transport or Earaday rotation experiment 
on films. Conceivably, the AC-field-driven AHE may ex¬ 
ist in other materials possessing non-trivial band topol¬ 
ogy such as ferromagnetic metals jsij . 
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Supplemental Material: 

When Chiral Photons Meet Chiral Fermions— 
Photoinduced Anomalous Hall Effects in Weyl Semimetals 

Ching-Kit Chan, Patrick A. Lee, Kenneth S. Burch, Jung Boon Han, and Ying Ran 


EFFECTIVE WEYL-FLOQUET HAMILTONIAN 

Here, we provide more details about the derivation of the effective Weyl-Floquet Hamiltonian in the main text. We 
start with a general Y-band tight-binding model H{k) of a Weyl semimetal (SM). For convenience, we choose the 
energy eigenstates at kw as the basis and shift the energy so that the Weyl node is at zero energy. The linearized 
Hamiltonian iJun in this basis has the follow block form: 

U^H{q + kw)U = Hnn + 0{q^) = 

where the columns of U are the eigenstates of H{kw) with the first two columns corresponding to the two Weyl 
bands. We introduce a real 3-vector a and a 3 x 3 real matrix [3 to parameterize the Weyl band structure. Here, the 
Hermitian matrices Di correspond to the {N — 2) high energy bands and Dq is diagonal, while Ci are 2 x (Y — 2) 
matrices describing the mixing between the Weyl and high energy bands. 

Consider a electromagnetic field propagating along the ^-direction. After coupling H\i^ with the gauge field: 
q ^ q-\- eA(t)^ with A{t) = {A^ cos(ci;t). Ay sin(ci;t -h 0), 0), we obtain: 

Hun ^ i^iin + ( 18 ) 


with H_i = 


and 


n = - f~ + Ax(3xj(7j - ie^^AyPyjdj A^C^ - ie'^^AyCy\ 

' “ 2 V ^xCl - ie^^AyCl A,D, - ie^^AyDy) 

In the perturbative regime A/uj 1, the time-independent effective Hamiltonian near kw is given by [l[: 

UJ 

_( 0{q)+0{A'^) \ 

[0{q)+0{A^) Do + 0{q) + 0{A^))’ 


(19) 


( 20 ) 


with 


Hwf{^ =qiaiao + qiPi^aj - [tijkP.iPyjak - Uc,Cl - CyCD], ( 21 ) 

which is the result in the main text. Perturbation theory tells that Ywf(^) is the correct 2-band effective Weyl- 
Floquet Hamiltonian up to 0{q‘^) -\-0{qA‘^) -f-O(A^). Note that the [3 term inside the bracket can be understood from 
a 2-band itself, while the Ci terms are contributed by the linear-g' mixing between the Weyl and high energy 

bands. 


CALCULATIONS FOR THE DRIVEN LATTICE MODEL 

In this section, we detail the calculation for the Weyl model studied in the main text. We adopt the four-band 
model proposed by Ojanen [H on a diamond lattice with an inversion symmetry breaking field. In momentum space, 
the undriven Hamiltonian is given by: 

5 

iJo = Y - *er 4 r 5 , 

i=l 


( 22 ) 
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where are Gamma matrices acting on the spin and orbital degrees of freedom. We take Fi = a F 2 = 
F 3 = cFz ^ Tz^ F 4 = (Tx and F 5 = (7y. Here, 


Di (k) = 2 A sin — cos-cos — 


kz 


ko. 




D 2 {k) = 2 A sin cos —— cos 


kcr. 


k. 


Ds{k) = 2 A sin — ( cos — — cos 


D 4 {k) = t ( 1 + cos 


“ 1 “ 


■ COS ■ 


kx H“ kz 


cos 


ky + kz 


D^{k) = t ( sin 


. kx “ 1 “ ky . kx “h kz . ky “h k^ 


+ sm - 


+ sm ■ 


(23) 


where t, A, and e parameterize the nearest-neighbor hopping, spin-orbit coupling and the inversion breaking potential, 
respectively. When |e| < 4|A|, this system possesses 6 pairs of Weyl nodes [ 2 | at kw = (i ^05 05 27r), (±/co, 27r, 0), 
(0, zb/co, 27r), (27r, ±/co, 0), (27r,0,±/co) and (0, 27r, ±/co), with ko = 2 sin“^[e/(4A)] (see Fig. 2(a) in the main text). 
However, due to the TR symmetry, the undriven system does not have a net anomalous Hall conductivity (i.e. 

^ = EiX^w^k\P = o). 

In the presence of a circularly polarized light, the Hamiltonian becomes time-dependent throng the Peierls sub¬ 
stitution k + H(cosco’t, sinu;t, 0). This allows the definition of the Floquet Hamiltonian [3|: iFi? |(/)^(t)) = 

[H{t) — idt] |0n(O) = where and |0n(^)) are the Floquet eigenenergies and eigenvectors satisfying the 

time-periodicity |0n(O) = |0n(^ + 27r/co’)). Fourier expanding |0n(^)) further enables us to express the Floquet Hamil¬ 
tonian as (n'|iFi?|n) = iFn-n' + nuj5n,n'h^ where \n) denotes the eigenmode with frequency nuu. 

The Fourier expansion of H{t) = are analytically tractable. The result is: 


5 

H-n = ^ di^-n{k)ri — ieF4F5(5n,0, (24) 

i=l 


where 

di-n{k) = A |2J„ ^ 


d2-n{k) = A - — Jn pT cos ^ sin 


cos — sm 
2 


A 


kx + nil 


d 3 -n{k) = 2\Jn — sin ^ 




2 

ky -b riTT 
2 

ky + nir 


-Jn 


( A 


in-K , f kx ^ ky ^ m:\ , . f kx — ky ^ nir 


\V2 


-b e 4 sm 


rvy 




_ iriTv f kx ky n7r\ iriTT f kx — ky niT 


e 4 sm 


e 4 sm 


Vx rvy 


kx + niT 


dA-n{k) = t\ 5nfi + Jn ( e ‘ 4 " cos - J 

kx + ky + riTT 


-HHL f kx + ky + mr 


Jn ( 4 


kx-\-kz-{- nir 


+ i-ir 


ky kz nir 


d5-n{k) =t\Jn(-^)e 4 sm 


7 / ^ 

+ Fn — 


sm 


kx + kz + nn 


+ {—i)^ sin 


ky kz + nir 


)]}■ 

(25) 


The Weyl node shifts in the main text (Fig. 2 and 3) are obtained by numerically solving the full Floquet Hamiltonian. 
For demonstration purposes, we used the parameter values t = l, e = 3, A = l,cj = l.l and A up to a/O^. Our 
findings are qualitatively unchanged in different parameter regimes, as long as the bare system possesses Weyl nodes 
(i.e. |e| < 4|A|) and the driving amplitude is not large enough to cause topological phase transitions by merging Weyl 
nodes. In our numerical calculations using A up to a/ 0~3, the Floquet spectrum converges rapidly by including up to 
n = ±2 Floquet eigenmodes. 


REALISTIC EXPERIMENTAL ESTIMATIONS 

We provide more supports and details for the estimated transport measurement discussed in the main text. First, let 
us define the spacing between the voltage contacts as current contacts as the sample thickness t and penetration 
depth 6. We assume the light equally illuminates the contact region with a power P. The sample has a refraction 
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index n{uj), a reflectance R{uj)^ a real part of the optical conductivity cri(ci;) and a dc longitudinal conductivity 
Furthermore we consider an experiment using a CO 2 CW laser {huo = 117 meV) and a constant current The 
measured Hall voltage is: 




y5/t 


X X It 


“t“ 1Ixt 


(26) 


where for simplicity we assume a uniform axy throughout the penetration depth, while the current is evenly spread 
through the entire thickness (hence the factor 6/t in axy). Since axx ^ we only need to be concerned with 
^xy^I{(^xx^)' Each Weyl node contributes a axy = note the relations ^ 

5{(jj) = Thus, accounting for the 24 Weyl nodes in TaAs, we find: 


T. 24 hvFe‘^ 2P(1 - i?) L ^ 

^ ^TT /fe N3 X TT — r\ X X 4 (27) 

a^^E 27rh {huj)^ lxly(^i{^) lx 

From previous measurements of the optical Q and DC conductivities [H and vp of TaAs [^-0, we take the values at 
300 K: P ^ 0.8, ai ^ 2 x 10^ axx ~ 1-7 x 10^ and an average hvp ~ 2 eVA. Considering a sample 

with t = 100 nm, lx = ly = 100 /im and using P=1 W, Ix=f A, the resultant Hall voltage is Vy ^ 130 nV, which is 
easily observable in realistic experiments. 

As we remark in the main text, the AHE is also measurable in a pump-probe experiment. Here, we give an order- 
of-magnitude estimation for axy in a pulsed laser experiment. Based on the Floquet topological insulator @, we take 
a pump laser energy huj ^ 0.12 eV and an electric field intensity eA ^ 0.01 — 0.02 together with the group 

velocity for TaAs used above, the resultant axy ^ 25 - 100 n-^cm-^ . 


DOPING EFFECTS 


Now, we explain the consequences of doping to the photoinduced AHE. Eor a better comparison, let us recap the 

findings for axy in the main text, where the chemical potential lies at the Weyl nodes (/i = 0). Before the drive, 

—*■ 2 

each Weyl node pair separated by K contributes a axy{lJi = A = 0) = ^^Kz. If the system has TR symmetry, the 
contributions from the TR Weyl node pairs cancel each other and the total axy{lJi = A = 0) = 0. When there is a 
circularly polarized drive, each Weyl node pair results in a photoinduced Saxy{fi = 0, A) ^ 2hvFe A being 
the Eermi velocity). The photoinduced contributions from different Weyl pairs always add. 

When the initial chemical potential is tuned away from the node, but still within the linearly dispersive regime, our 
result is unaffected. We first review the argument for the undriven case used in Ref. 0- In the linear regime, the 
low-energy Hamiltonian of a undriven Weyl node is i^eff = ' k)ai — //, where k is the momentum away from 

the node. The Berry curvature B{k) is odd about the node: 


B{-k) = -B{k), 

as can be shown by performing a formal TR operation: Ti7eff(^)T“^ = i7eff(—^)- As a result, we have 


^xy 


{li,A = 0) - (Ta;y(0, A = 0) 


A'w'^h 


I 


(fkB{k) ■ k^ 


0</c</eF 




(28) 


(29) 


where kp is the Eermi momentum and ci is some microscopic constant. In fact, an ideal Weyl node has B{k) (x k/k!^ 
and the above integral vanishes due to the cancellation between contributions from the -^kz and —kz planes. The 
0[/i^] correction comes from the departure from the ideal Weyl node. Note also that in a TR symmetric system, 
(^xyik'if^) is always zero. 

We now apply the same argument to the driven system. Since a driven Weyl node is still a monopole, Eq. (|28|) still 
holds and we have a similar relation: 

g2 

^xyik'') A) — axy{0^ A) ^ 2ph ^2 A )fi , (30) 

where the C 2 A^ term is a perturbative correction. Combining Eq. p?!i . (pn|i and Saxy{0, A), we have: 


daxyifJ^, A) = cr^^(/i. A) - cr^^(/i, 0) 


2ph \ 


f Tiwkvpe^ 2 

-h C2/i 


(31) 


UJ 
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FIG. 4. Demonstration of the doping effects of the driven two-band Weyl model (see. Eq. (l3^ b (a) Weyl-Floquet bands at 
kx — ky — 0. (b) axy(iJ^i,A) as the chemical potential fii is increased from 0 to 1.5. When fi is within the linearly dispersive 
regime, Saxy{/J^, A) oc A^. Parameters used: t = mi = m 2 = oj = 1. 


for systems with nw Weyl nodes. Therefore, as we tune /i away from the Weyl node, 6axy is still proportional to 
with a 0[jj!^] modification to the prefactor. 

When we go beyond the linearly dispersive region, contributions from large momentum enter and the AHE is 
further modified. Based on the discussion above, we would anticipate a nonlinear change: ^ /(//, 0) and 

C 2 A^(A A)|a=o^^ + •••, where / depends on microscopic band-structure details. 

To illustrate these effects, we numerically compute (Jxy{lJi^A) for a two-band model without TR symmetry 0: 


H{k) = 2t sin 


T 


(Jx + 2t sin 


kx ko^ 


(7y + [2mi — m 2 (cos kx + cos ky) 2 cos kz] < 


(32) 


We consider the parameters t = mi = m 2 = 1, at which H supports a Weyl node pair at /c = (0,0, ±7r/2), so that the 
undriven (Jxy{0A) = ^ 'tt. Figure HJ^a) and (b) present the Floquet bands and axy{fi,A), respectively, based on a 

perturbative calculation as // is tuned away from the node. The results show the expected behavior of 5(Txy{lJi^ A) oc A^ 
when {i is small (linearly dispersive regime), and nonlinear corrections for larger /i and A. 

We note that, in the calculation above, we have assumed only the bands below the chemical ^tential are filled. In 
a general driven system, excited states can also be populated and contribute to (Jxy (see e.g. [l|). While this could 
be an issue for the quantization of the Hall conductance, it does not qualitatively affect our results since we do not 
require a quantized conductance. The excited state population leads to a correction to the Berry curvature integral 
(like Eq. (|29|) ) and only modifies the proportionality between dcjxy and A^. 
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